Observation of Flux Reversal in a Symmetric Optical Thermal Ratchet 
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We demonstrate that a cycle of three holographic optical trapping patterns can implement a 
thermal ratchet for diffusing colloidal spheres, and that the ratchet-driven transport displays flux 
reversal as a function of the cycle frequency and the inter-trap separation. Unlike previously de- 
scribed ratchet models, the approach we describe involves three equivalent states, each of which 
is locally and globally spatially symmetric, with spatiotemporal symmetry being broken by the 
sequence of states. 



Brownian motion cannot create a steady flux in a sys- 
tem at equilibrium. Nor can local asymmetries in a static 
potential energy landscape rectify Brownian motion to 
induce a drift. A landscape that varies in time, however, 
can eke a flux out of random fluctuations by breaking 
spatiotemporal symmetry Q, 0, 0, • Such flux-inducing 
time-dependent potentials are known as thermal ratch- 
ets H, E|, an d their ability to bias diffusion by rectify- 
ing thermal fluctuations has been proposed as a possible 
mechanism for transport by molecular motors and is be- 
ing actively exploited for macromolecular sorting • 

Most thermal ratchet models are based on spatially 
asymmetric potentials. Their time variation involves dis- 
placing or tilting them relative to the laboratory frame, 
modulating their amplitude, changing their periodicity, 
or some combination, usually in a two-state cycle. Chen 
demonstrated that a spatially symmetric potential still 
can induce drift, provided that it is applied in a three- 
state cycle, one of which allows for free diffusion |S| . This 
idea since has been refined $] and generalized . 

The space-filling potential energy landscapes required 
for most such models pose technical challenges. Further- 
more, their relationship to the operation of natural ther- 
mal ratchets has proved difficult to establish. 

This Letter describes the first experimental demon- 
stration of a spatially symmetric thermal ratchet, which 
we have implemented with holographic optical traps 
[ill H2I fl3| . The potential energy landscape in this sys- 
tem consists of a large number of discrete optical tweezers 
[l3 | , each of which acts as a symmetric potential energy 
well for nanometer- to micrometer-scale objects such as 
colloidal spheres. We arrange these wells so that colloidal 
spheres can diffuse freely in the interstitial spaces but are 
localized rapidly once they encounter a trap. A three- 
state thermal ratchet then requires only displaced copies 
of a single two-dimensional trapping pattern. Despite its 
simplicity, this ratchet model displays flux reversal 0,0] 
in which the direction of motion is controlled by a bal- 
ance between the rate at which particles diffuse across 
the landscape and the ratchet's cycling rate. 

Often predicted, and inferred from the behavior of 
some natural molecular motors and semiconductor de- 
vices , flux reversal has been directly observed in corn- 
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FIG. 1: A spatially-symmetric three-state ratchet potential 
comprised of discrete potential wells. 

paratively few systems. Previous demonstrations have 
focused on ratcheting of magnetic flux quanta through 
type-II superconductors in both the quantum mechanical 
1161 and classical regimes, or else have exploited the 
crossover from quantum mechanical to classical trans- 
port in a quantum dot array |l8j . Unlike the present 
implementation, these exploit spatially asymmetric po- 
tentials and take the form of rocking ratchets ||. A 
massively parallel hydrodynamic ratchet for fluid-borne 
objects driven by oscillatory flows through asymmetric 
pores also shows size-dependent flux reversal Eflj . 
In this case, however, the force field is provided by the 
divergence-free flow of an incompressible fluid and so 
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cannot be described as arising from a potential energy 
landscape. Rather, this is an instance of a so-called 
drift ratchet |20| . Other pioneering implementations of 
classical force-free thermal ratchets also were based on 
asymmetric potentials, but did not exhibit flux reversal 

hi mm mm. 

Figure^shows the principle upon which the three-state 
optical thermal ratchet operates. The process starts out 
with a pattern of discrete optical traps, each of which 
can localize an object. The pattern in the initial state 
is schematically represented as three discrete potential 
energy wells, each of width a and depth Vo, separated by 
distance L. A practical trapping pattern can include a 
great many optical traps organized into manifolds. The 
first pattern is extinguished after time T and replaced 
immediately with the second, which is displaced from 
the first by L/3. This is repeated in the third state with 
an additional step of L/3, and again when the cycle is 
completed by returning to the first state 

If the traps in a given state overlap those in the state 
before, a trapped particle is transported deterministically 
forward. Running through this cycle repeatedly transfers 
the object in a direction determined unambiguously by 
the sequence of states, and is known as optical peristalsis 
|2(j|. The direction of motion can be reversed only by 
reversing the sequence. 

The optical thermal ratchet differs from this in that 
the inter-trap separation L is substantially larger than 
a. Consequently, particles trapped in the first pattern 
are released into a force-free region and can diffuse freely 
when that pattern is replaced by the second. Those par- 
ticles, such as the example labeled "forward" in Fig. ^ 
that diffuse far enough to reach the nearest traps in the 
second pattern rapidly become localized. A comparable 
proportion of this localized fraction then can be trans- 
ferred forward again once the third pattern is projected, 
and again when the cycle returns to the first state. 

Unlike optical peristalsis, in which all particles are pro- 
moted in each cycle, the stochastic ratchet transfers only 
a fraction. This, however, leads to a new opportunity. 
Other particles that are too slow to catch the forward- 
going wave might still reach a trap on the opposite side 
of their starting point while the third pattern is illumi- 
nated. These particles would be transferred backward by 
L/3 after time 2T, as shown in the "reverse" trajectory 
in Fig.Hl 

For particles of diffusivity D, the time required to 
diffuse the inter-trap separation is r = L 2 /(2D). If 
we assume that particles begin each cycle well localized 
at a trap, and that the traps are well separated com- 
pared to their widths, then the probability for ratch- 
eting forward by L/3 during the interval T is roughly 
P F w exp(-(L/3) 2 /(2L>T)), while the probability of 
ratcheting backwards in time 2T is roughly Pr w 
exp(— (L/3) 2 /(4L>T)). The associated fluxes of particles 
then are v F = P F L/(3T) and v R = -PrL/(6T), with the 



dominant term determining the overall direction of mo- 
tion. Crudely, then, we expect the direction of induced 
motion to reverse when T/r < (18hi2) _1 w 0.08. 

More formally, we can model an array of evenly spaced 
optical traps in the u-th pattern as Gaussian potential 
wells 
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where n = 0, 1, or 2, and N sets the extent of the land- 
scape. The probability density p(x,t)dx for finding a 
Brownian particle within dx of position x at time t in 
state n evolves according to the master equation 

p(y, t + T) = J P n (y, T\x, 0) p(x, t) dx, (2) 

characterized by the propagator 

P n {y,T\x,Q) = e L ^ T 5{y-x), (3) 
where the Liouville operator for state n is 

Ln(y) = D (J^ - pl-V^y)) , (4) 

with V^(y) = and where /3 _1 is the thermal energy 
scale. 

The master equation for a three-state cycle is 

p(y, t + 3T) = J P 123 (y, 3T\x, 0) p(x, t) dx, (5) 
with the three-state propagator 

J Pi 23 (2/,3T| a; ,0) = J d Vl dy 2 P 3 (y,T\y 2 ,0)x 

P 2 (y 2 ,T\y 1 ,0)P 1 (y 1 ,T\x,0). (6) 

Because the landscape is periodic and analytic, Eq. (J5J 
has a steady-state solution such that 



p(x, t + 3T) = p(x, t) 



(7) 
(8) 



The mean velocity of this steady-state then is given by 



Pi 23 {y,3T\x,0) 
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3T 



pi 23 (x)dxdy, (9) 



where Pi 23 {y, 3T\x, 0) is the probability for a particle 
originally at position x to "jump" to position y by the end 
of one compete cycle, (y — x)/ (3T) is the velocity associ- 
ated with making such a jump, and pi 23 (x) is the fraction 
of the available particles actually at x at the beginning 
of the cycle in steady-state. This formulation is invariant 
with respect to cyclic permutations of the states, so that 
the same flux of particles would be measured at the end 
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of each state. The average velocity v therefore describes 
the time-averaged flux of particles driven by the ratchet. 

Figure^a) shows numerical solutions of this system of 
equations for representative values of the relative inter- 
well separation L/a. If the interval T between states 
is very short, particles are unable to keep up with the 
evolving potential energy landscape, and so never travel 
far from their initial positions; the mean velocity van- 
ishes in this limit. The transport speed v also vanishes 
as 1/T for large values of T because the induced drift 
becomes limited by the delay between states. If traps 
in consecutive patterns are close enough (L = 6.5 a in 
Fig. Eta)) particles jump forward at each transition with 
high probability, yielding a uniformly positive drift ve- 
locity. This transfer reaches its maximum efficiency for 
moderate cycle times, T/r w 2V2{L/a)(/3V )~ 1 . More 
widely separated traps (L = 13 a in Fig.^a)) yield more 
interesting behavior. Here, particles are able to keep up 
with the forward-going wave for large values of T. Faster 
cycling, however, leads to flux reversal, characterized by 
negative values of v. 

We implemented this thermal ratchet protocol for 
a sample of 1.53 /zm diameter colloidal silica spheres 
(Bangs Laboratories, lot number 5328) dispersed in wa- 
ter, using potential energy landscapes created from ar- 
rays of holographic optical traps |llj, [l2j, |l3J ■ The sam- 
ple was enclosed in a hermetically sealed glass cham- 
ber roughly 40 /im thick created by bonding the edges 
of a coverslip to a microscope slide and was allowed 
to equilibrate to room temperature (21 ± 1°C) on the 
stage of a Zeiss S100TV Axiovert inverted optical mi- 
croscope. A 100 x NA 1.4 oil immersion SPlan Apo ob- 
jective lens was used to focus the optical tweezer array 
into the sample and to image the spheres, whose mo- 
tions were captured with an NEC TI 324A low noise 
monochrome CCD camera. The micrograph in Fig. Hfb) 
shows the focused light from a 20 x 5 array of opti- 
cal traps formed by a phase hologram projected with 
a Hamamatsu X7550 spatial light modulator [2s| . The 
tweezers are arranged in twenty-trap manifolds 25 /im 
long separated by Lq — 6.7 /im. Each trap is pow- 
ered by an estimated 2.5 ± 0.4 mW of laser light at 
532 ran. The particles, which appear in the bright-field 
micrograph in Fig.^c), are twice as dense as water and 
sediment to the lower glass surface, where they diffuse 
freely in the plane with a measured diffusion coefficient 
of D = 0.33 ± 0.03 /im 2 /sec, which reflects the influence 
of the nearby wall. Out-of-plane fluctuations were mini- 
mized by projecting the traps at the spheres' equilibrium 
height above the wall |3JJ] . 

We projected three-state cycles of optical trapping pat- 
terns in which the manifolds in Fig. Hfb) were displaced 
horizontally by — Lq/3, 0, and Lq/3, with inter-state de- 
lay times T ranging from 0.8 sec to 10 sec. The parti- 
cles' motions were recorded as uncompressed digital video 
streams for analysis [2^. Between 40 and 60 particles 
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FIG. 2: (a) Crossover from deterministic optical peristalsis at 
L = 6.5 a to thermal ratchet behavior with flux reversal at 
L — 13 a for a three-state cycle of Gaussian well potentials at 
/3Vo = 8.5, a = 0.53 ^m and D — 0.33 /im 2 /sec. Intermediate 
curves are calculated for evenly spaced values of L. (b) Image 
of 20 x 5 array of holographic optical traps at Lo = 6.7 yxa. 
(c) Image of colloidal silica spheres 1.53 ^m in diameter in- 
teracting with the array, (d) Rate dependence of the induced 
drift velocity for fixed inter-trap separation, Lq- (e) Separa- 
tion dependence for fixed inter-state delay, T = 2 sec. 

were in the trapping pattern during a typical run. so that 
roughly 40 cycles sufficed to acquire reasonable statistics 
under each set of conditions without complications due 
to collisions. We also tracked particles outside the trap- 
ping pattern to monitor their diffusion coefficients and 
to ensure the absence of drifts in the supporting fluid. 
The results plotted in Fig. E2d) reveal flux reversal at 
T/t w 0.03. Excellent agreement with Eq. © is ob- 
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tained for /3V = 8.5 ± 0.8 and a = 0.53 ± 0.01 (im. 

The appearance of flux reversal as one parameter is 
varied implies that other parameters also should control 
the direction of motion [5|. Indeed, flux reversal is ob- 
tained in Fig. Hje) as the inter-trap separation is var- 
ied from L = 5.1 /im to 8.3 /jm at fixed delay time, 
T = 2 sec. These results also agree well with predictions 
of Eq. (PJ, with no adjustable parameters. The same 
effect also should arise for different populations in a het- 
erogeneous sample with different values of D, Vo and a 
|3l|,|32j. In this case, distinct fractions can be induced to 
move simultaneously in opposite directions. 

Such sensitivity of the transport direction to details 
of the dynamics also might play a role in the function- 
ing of molecular motors such as myosin- VI whose ret- 
rograde motion on actin filaments compared with other 
myosins has excited much interest l33|. This molecular 
motor is known to be nonprocessive 34] ; its motion in- 
volves a diffusive search of the actin filament's potential 
energy landscape, which nevertheless results in unidirec- 
tional hand-over-hand transport . These characteris- 
tics are consistent with the present model's timing-based 
flux reversal mechanism, and could provide a basis to 
explain how small structural differences among myosins 
could lead to oppositely directed transport. 
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